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ABSTRACT 

> , 

The oscillation spectrum of v Eri is the richest known for any variable of the 
■ [3 Cephei type. We interpret the spectrum in terms of normal mode excitation and 

construct seismic models of the star. The frequency data combined with data on mean 
colours sets the upper limit on the extent of overshooting from the convective core. 
We use data on rotational splitting of two dipole {1 = 1) modes (gi and pi) to infer 
\ properties of the internal rotation rate. Adopting a plausible hypothesis of nearly 

— . ■ uniform rotation in the envelope and increasing rotation rate in the /i-gradicnt zone, 

*-£h ' we find that the mean rotation rate in this zone is about three times faster than in the 

envelope. In our standard model only the modes in the middle part of the oscillation 
Q ■ spectrum are unstable. To account for excitation of a possible high-order g-mode at 

S-^ | v = 0.43 cd 1 as well as p-modes at v > 6 cd _1 we have to invoke an overabundance 

r/3 , of Fe in the driving zone. 
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1 INTRODUCTION seismic models of a star, that is, the models constrained by 

the frequency data. 

In this work, we attempt to make the best use of the 
frequency data to address unsolved problems in stellar evolu- 
The rich oscillation spectrum of v Eri obtained in a re- tion theor y concerning element mixing in convectively stable 
cent multisite campaign (Handler et al. 2004, Aerts et al. la y ers and angular momentum evolution. The problems are 
2004) holds the best prospects for seismic sounding of the related because rotation induces fluid flows that may cause 
interior of a B-star but also presents a considerable chal- mixing. It is important to disentangle effects of overshooting 
lenge to stellar pulsation theory. The sounding is facilitated tom the convective core and effects of rotationally induced 
by the fact that for several excited modes the spherical har- element mixing. The first effect depends only on stellar mass, 
monic indices were determined (De Ridder et al. 2004). The The second one depends on randomly distributed stellar an- 
challenge is to explain mode excitation in an unexpectedly S ular momenta. In this context, v Eridani, which is a very 
broad frequency range slowly rotating B star, yields a useful extreme case. How- 

ever, it would be important to try to obtain similar quality 
data for more rapidly rotating objects. 

We are very curious how the angular velocity of rota- 
The oscillation spectrum of v Eri is shown in Fig. 1. Be- tion behaves in stellar interiors. Let us remind that, while 
fore the campaign only the I = mode and the (i = 1, gj) the results of helioseismic sounding of radial structure essen- 
triplet were known. However, for the triplet, the identifica- tially confirmed the standard solar model, the results for the 
tion of the spherical harmonic was uncertain. Now, thanks to internal rotation rate came as a surprise. The outstanding 
the analysis of accurate multicolour photometry and high- 
resolution spectroscopy, the uncertainty was eliminated. Re- 
liable information about i is a crucial input for constructing * E-mail: alosza@camk.edu.pl 
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Figure 1. Oscillation spectrum of v Eri. The peaks repre- 
sented with the broken vertical lines are regarded uncertain. 
The four modes known before the campaign are between v = 
5.6 and 5.8 cd" 1 . Dc Ridder et al. (2004) identified them as an 
t = 1 triplet and an £ = singlet. The newly discovered triplet 
at v « 6.25 cd" 1 is identified as £ = 1 (De Ridder et al. 2004). In 
this paper we will show that the radial mode is the fundamental 
and the two dipole modes are gi and pi , respectively. Only these 
mode frequencies are used in our seismic sounding. 



question in the case of massive stars is whether a substantial 
radial gradient of the angular rotation rate may be present 
in their interiors. The question is related to another very 
interesting and important problem which is the origin of 
magnetic fields in B type stars. 

A seismic model of a pulsating star should not only ac- 
count for the measured frequencies but also for the instabil- 
ity of the detected modes. The latter requirement sets a dif- 
ferent kind of constraint than the former as it concerns only 
properties of the outer layers, where most of contribution 
to mode driving and damping arises. Since the frequency 
spectrum of v Eri seems unusually broad for an object of its 
type, constraints in this case are particularly interesting. 

The structure of this paper is as follows. In Sect. 2 
we construct models of the internal structure for v Eri 
and calculate its oscillation frequencies. We find models for 
which the oscillation frequencies of three m = modes - 
(£ = 0, p x ), (£ = 1, g]) and (£ = 1, p x ) - reproduce the 
observations as well as have effective temperatures and lu- 
minosities within the observational error box. In Sect. 3 we 
use these models and data on the triplet structure to infer 
the internal rotation of the star. 

The problem of how to make the observed modes pul- 
sationally unstable is discussed in Sect. 4. The tasks of con- 
structing the model of v Eri whose mode frequencies match 
the data and of the identification of the driving effect are not 
quite independent. The driving effect in B stars strongly de- 
pends on the iron content and its distribution in outer layers. 
In Sect. 5, we study how modification of the iron distribution 
affects the frequencies. We also propose plausible identifica- 
tions of all modes detected in the v Eri oscillation spectrum. 

2 MODEL OF THE STELLAR INTERIOR 

For the purpose of constructing models of v Eri we ignored 
all effects of rotation which, as we have checked, was fully 



justified. Accordingly, we make use here only of the m = 
mode frequencies despite, in general, rotation causes a shift 
of the centroid frequency. We consider models of massive 
main sequence stars in the phase of radius expansion. A 
model is characterised by its mass, M, effective tempera- 
ture, Tog, initial hydrogen, X, and heavy element, Z, abun- 
dance, as well as the overshooting parameter, a ov . Here we 
assume the standard solar mixture of heavy elements (de- 
partures will be considered in section 5) and fix X at 0.7. 
For opacity and the equation of state we use the OPAL tab- 
ular data (Iglesias & Rogers 1996, Rogers et al. 1996). The 
nuclear reaction rates are the same as used by Bahcall and 
Pinsonneault (1995). 



2.1 v Eridani's position in the theoretical H-R 
diagram 

The position of the star in the log T e g — log L diagram is not 
very accurately known. Our temperature determination is 
based on Stromgren photometry and on the tabular data by 
Kurucz (1998). For L we use the Hipparcos parallax, Lutz- 
Kelker correction to the absolute magnitude (Lutz & Kelker 
1973), and the bolometric correction from Kurucz's data. In 
this way we obtained the following numbers 

log T cB = 4.346 ±0.011 and log L = 3.94 ± 0.15. 

For log T e ff we include measurement errors and uncertainty 
of calibration. For logL we take into account errors of the 
parallax and uncertainty of the bolometric correction. Our 
numbers are similar to those quoted by Dziembowski & 
Jerzykiewicz (2003) and will be updated by De Ridder et 
al. (2004). 



2.2 Identification of the radial overtone of the 
observed modes 

Before starting seismic modelling of v Eri, we must complete 
the mode identification. The work by De Ridder et al. (2004) 
provided unambiguous I identifications for seven modes, and 
m determinations by inference, but the radial overtone of 
the modes has not been determined. Fortunately, one of the 
detected pulsation modes is radial, which eases this task 
considerably. 

Consequently, we follow Dziembowski & Jerzykiewicz 
(1996) and only consider two one-dimensional families of 
models. The position of v Eri in the HR diagram is only con- 
sistent with the radial mode at 5.7634 cd _1 being the funda- 
mental or the first overtone (as indicated later in Fig. 3) . For 
these two possibilities, we examine the mass dependence of 
theoretical I = 1 mode frequencies of standard models and 
match it to the observed I = 1 mode frequencies (Fig. 2) . 

Because of the mixed-type character of some of the 
modes we can only match the observed frequencies of the 
two lowest-order 1=1 modes within certain mass ranges. 
However, under the assumption of the radial mode being the 
first overtone, no acceptable fit can be found for the I = 1 
mode of highest frequency. Hence we can eliminate this pos- 
sible identification (which would also be less consistent with 
the star's position in the HR diagram) and identify the radial 
pulsation mode of v Eri as the fundamental. Consequently, 
the £ = 1 modes are the first g-mode (gi) and the first two 
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Table 1. Parameters of fitted models. The symbols denote : X c - hydrogen abundance in the core, M mlx - mass of the mixed core, 
which is the same as mass of the convective core at a ov = 0, r c - radius of the convective core, r c o - radius of the convective core at 
ZAMS, which corresponds to the top of the /^-gradient zone. 
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Figure 2. Identification of the radial overtone of the observed 
modes of v Eri. We assume the radial mode to be the fundamen- 
tal in the upper panel and to be the first overtone in the lower 
panel. The mass dependence of theoretical t = 1 mode frequencies 
(dashed lines) is compared to the observed 1 = 1 mode frequen- 
cies (horizontal lines). We find two ranges in mass (denoted by 
the vertical dotted lines) where the two 1 = 1 modes of lowest fre- 
quency match the observed frequencies well. However, the 1 = 1 
mode near 7.9 cd _1 can only be reproduced in the case of the 
radial mode being the fundamental. 



p-modes (pi and P2). We note that the theoretical fit of the 
standard model to the observed p 2 mode is also not very 
satisfactory at first glance; we will return to this problem 
and its possible solution later. 



2.3 Fitting frequencies of the I = and two 
I = l,m = modes 

We calculated frequencies of these three modes in evolu- 
tionary sequences of stellar models characterised by the pa- 
rameters M, Z, a ov . The first two of them and T e g were 
regarded as adjustable parameters to fit the measured fre- 
quencies. The adopted precision of the fit was 10" 3 . The 
precision of the data is much higher but matching to such 
a precision is neither interesting nor reasonable in view of 
the uncertainties of the model calculations, such as the lim- 
ited accuracy of the opacity data. The upper limit of the 
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Figure 3. The evolutionary tracks for the models fitting frequen- 
cies of the (I = 0, pj), (i = 1, gj) and (ft = 1, p x ) modes in v 
Eri. Models were calculated with chemical composition param- 
eters Z = 0.015 (adjusted) and X = 0.7 (fixed). The values of 
M and a ov arc indicated. The values of T e fj and log L inferred 
from photometry and the Hipparcos parallax are shown with the 
error bars. Two dotted lines connect models with 5.7633 c/d as 
the radial fundamental (pi) and first overtone (p 2 ) mode. 



a ov parameter was derived from the bounds on the effective 
temperature. 

We obtained models with a ov = 0.0 and q ov = 0.1 
which satisfy all the observational constraints. The inferred 
value of the metallicity parameter is Z = 0.0150±0.0001, the 
same at a ov = 0.0 and 0.1. Parameters of the two selected 
models are given in Table 1. 

Both models are well within the error box in the the- 
oretical H-R diagram shown in Fig. 3. Increasing a ov above 
0.12 would place the seismic model below the lower bound 
of T e ff. We, thus, conclude that the data on v Eri are consis- 
tent with negligible overshooting distance and set an upper 
bound on it at ev ov = 0.12. Further constraining the extent of 
overshooting may be possible when the effective temperature 
determination is improved. Our constraint is in agreement 
with a ov = 0.10 ± 0.05 derived by Aerts et al. (2003) from 
frequency analysis for the (3 Cephei star HD 129929, which 
is also slowly rotating. 

The frequency distance between the (I = 1, g x ) and 
(£ = 0, pj) modes is indeed a very sensitive probe of over- 
shooting. This is illustrated in Fig. 4, where we show the 
evolution of the three mode frequencies as a function of T e g . 
This is the main application of the (£ = 1, g x ) mode. How- 
ever, without the (£ = 1, p x ) mode frequency the assessment 
of a ov would not be possible unless we have an accurate ob- 
servationally determined Z value which we do not. 



4 A. A. Pamyatnykh et al. 



6.5 



5.5 



V s 

X 


— ^< 




i i | i i_ 
a ov =0.0 






<^ 














' — s 




iil,, - 



5.5 



: v, [c/d] 


X 1 V 

X? X 

V s - 

' ^ 




1 [ I |_ 

a ov =0.1 






4"^ 





















4.37 



4.36 



4.35 
log T. 



4.34 



4.33 



Figure 4. The evolution of the calculated frequencies of the se- 
lected modes in the model sequences shown in Fig. 2. Two dotted 
vertical lines mark the allowed T e fj range. Three horizontal lines 
correspond to the measured frequencies. Solid vertical lines cor- 
respond to T e ff for the fitted models. The fittings are marked by 
open circles. Note that the gi mode frequency comes close to the 
radial mode frequency in a very narrow range of T c ff and the lo- 
cation of this range is very sensitive to o v Crossing of the two 
frequencies takes place at the point of the minimum frequency 
distance between the two I = 1 modes (avoided crossing). Before 
the avoided crossing, the (I = 1, Pi) mode is nearly a pure acous- 
tic mode with its frequency evolution line being nearly parallel to 
the corresponding line for the radial mode. The ratio of the two 
frequencies depends on Z. 



2.4 Probing properties of low order I = 1 modes 

It is remarkable how much the three mode frequencies tell us 
about the star's internal structure. Identification of a radial 
mode is important because its frequency yields an accurate 
constraint on the model parameters; however, much more 
interesting is the information contained in the I = 1 modes. 

The probing property of a mode depends on the dis- 
tribution of its oscillation energy, E, in the stellar interior. 
The plots in the upper panel of Fig. 5 show £, the energy 
derivative with respect to the fractional radius, for the four 
consecutive 1=1 modes. In spite of rather close frequen- 
cies, modes gi and pi have a grossly different distribution of 
energy. This explains the differences in their evolution seen 
in Fig. 4 as well as the difference in the probing properties. 
One should also notice the peculiar £ (r) dependence for the 
gi and pi modes in the deep interior. A maximum (abso- 
lute in the former and local in the latter case) is reached 
at the boundary of the convective core where the Brunt- 
Vaisala, frequency has a derivative discontinuity. The g2 and 
all higher order g-modes have the maxima within the prop- 
agation zone. For our models it is the gi mode whose fre- 
quency is most sensitive to overshooting. In more evolved 
models this property will shift to pi . Dziembowski & Pamy- 
atnykh (1991) observed this property of certain low order 
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Figure 5. Upper panel: Oscillation energy distribution, £, for low 
order 1 = 1 modes in the v Eri model calculated with ct ov = 0. The 
plots are very similar in the models with a ov = 0.1. Lower panel: 
The Brunt- Vaisala frequency, i^bVi an d the Lamb frequency at 
t = 1, j^li, in the same model. 



g- modes, called by them g c , in 5 Scuti star models, and em- 
phasized the importance of detecting such a mode in the 
context of the overshooting problem, v Eri is the first object 
in which a g c mode was definitely detected. More than that, 
its rotational splitting has been measured. 

v Eridani is not exceptional amongst /3 Cep stars, which 
are 9 — 12M Q stars in the advanced main sequence phase of 
evolution. In models of such stars we typically find that the 
(I = 1, g x ) and (I = 1, pj) modes are both unstable. There 
are chances that they will be detected in a number of other 
stars of this type. The good luck in the case of v Eri is that 
we capture it before the avoided crossing of the i = 1 modes. 
After it, the data on the modes yield much less orthogonal 
information. 



3 INTERNAL ROTATION 
3.1 The I = 1 triplets 

The frequency structure of a rotationally split 1=1 mode 
may be characterised by the mean splitting S = 0.5(^+i — 
V-i) and the asymmetry A = v-\ + v+\ — 2vq. If all effects 
of higher order in the angular rotation rate, Q, as well as 
effects of a magnetic field are negligible, then A = 0. Both 
parameters, A and S, for the gi mode have been determined 
before the present campaign (Dziembowski & Jerzykiewicz 
2003). 

If we allow only an r-dependence in f2 then 



vo + 



171 f Rdr rr> , n , 2 n 
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The rotational splitting kernel, fC, is given by 
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Figure 6. The rotational splitting kernel for the gi and pi modes 
in the model calculated with o ov = 0. The vertical line at r c o/R 
marks the top of the /i-gradient zone. Note the difference between 
K, and E, shown in Fig. 5. Rotation within the convective core has 
hardly any effect on the i = 1 splitting. 
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where £ r and £/, are defined by the following expression for 
the displacement eigenvector 

c = (^e r + av h )yr. 

The kernels for the two modes considered are shown in Fig. 6. 
The quantities Do and D\ are quadratic in Q. The asymme- 
try of the triplet arising from the quadratic effect of rota- 
tion is 2D\. The quantity Do represents the frequency shift 
of the m — mode. As we have mentioned in section 2.1, 
in principle, this shift should be subtracted from the mea- 
sured frequencies of m = modes when we fit frequencies 
calculated for a spherical model. 

From the frequencies listed by Handler et al. (2004) for 
the gi mode we derive S g = 1.69 x 10~ 2 cd -1 (from here 
on we use subscript g to denote quantities referring to the 
gi mode, those referring to the pi mode will be subscribed 
with p), which is nearly the same as before (Dziembowski 
& Jerzykiewicz 2003). There is, however, a significant dif- 
ference in the asymmetry for which the previous value was 



A a 



-7.1 x 10 cd , but from the most recent data we 



now determine a value which is by more than a factor of 
two smaller. The former value was larger than calculated 
assuming uniform rotation with the rate inferred from S g . 
Dziembowski & Jerzykiewicz (2003) considered two possi- 
ble solutions. One was faster rotation in the envelope than 
in the interior, the second was a contribution to D\ from 
a hypothetical magnetic field. The current data taken in a 
158 day-long campaign and their time distribution are in- 
sufficient for a precise determination of A g . Thus we will not 
use this parameter here and we will not consider effects of a 
possible magnetic field. We will postpone interpretation of 
the A g value until we have more precise determination from 
the future data analysis. 

For the pi mode we have reliable frequencies only for the 
m — and +1 components. The frequency of the m = — 1 



component may be estimated from the value of A g , assuming 
that D\ arises from the second order effect of rotation with 
the same rate for both modes. In this way we get (assuming 
A g from earlier data which we regard as more reliable) 



Ap — Ag 



Di, P 
Di, g 



-1.13 x lO^ cd" 



and 



= 2i/ , p — v +i,p + A p — 6.2250 cd 



A g — will make very small difference. In a similar man- 
ner we find a frequency shift D , P = 7 x 10~ 4 for the m = 
mode. The number is below the adopted precision of the 
frequency fit in our model construction. For the linear rota- 
tional frequency splitting we get 

S p = v + i, p -v ,p-0.5A p = 0.01797 + 0.00056 = 0.0185cd _1 



3.2 Rotation rate in the envelope and in the 
/i-gradient zone 

Since we have data on rotational splitting only for two 
modes, our inference on the internal rotation must rely on 
simplifying assumptions. Results from seismic sounding of 
the solar internal rotation suggest that within chemically ho- 
mogeneous radiative layers the rotation rate should be close 
to uniform. A possible steep gradient may occur only in the 
chemically inhomogeneous zone around the convective core, 
where the /x-gradient stabilises differential rotation. With 
this in mind, we write 



and 



where 



QirSg — K c ,g^c ~f" Ke, g Qe 



K c,j = / K& >i = / 

JO J r c o 

and j = p or g. In our fitted models we have 

K c>g = 0.187(0.168), K e , g = 0.345(0.395) 

K c ,p = 0.092(0.101), K e> p = 0.718(0.682), 

where numbers in brackets refer to the model calculated with 
a ov = 0.1. The results are nearly the same for both models. 
We find 



3Q e . 



Figure 6 shows the averaging kernels for fl The fact 
that the convective core virtually does not contribute to the 
splitting is an exclusive property of £ = 1 modes which may 
be easily seen by considering the behaviour of eigenfunctions 
near the r — > singularity. For finite £'s we have then £ r — ► 
££h, implying (see Eq.(l)) K. — > if I — 1. Thus, the value 
Q c refers only to the /i gradient zone. It is clear that our 
result implies a very sharp decline of the rotation in the 
layer between the current core boundary and that at the 
ZAMS phase. How sharp the decline is depends on the form 
of the Q(r) dependence. A factor ~ 5 is derived if we assume 
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a linear decrease of D. between r = r c and r c o and constant 
Q for r > r c0 . 

With this simple form of fi(r) we find an equatorial 
velocity of rotation v e = 6.1 km/s, which is 2/3 of the 
value inferred from S g on the assumption of uniform ro- 
tation. Since the quadratic effect of rotation is still expected 
to arise mainly above r — r c , we have a problem with ac- 
counting for the pre-campaign value of A g . The values of 
the asymmetries calculated ignoring the contribution from 
below r = r c o are A g = —2.0 x 1CF 4 and A p = —2.6 x 1CT 4 . 
The time for a new discussion of the triplet asymmetry will 
come when we have more precise frequencies of the modes 
in both triplets. 



4 HOW THE MODES EXCITED IN v ERI ARE 
DRIVEN 

Oscillation driving in B-type stars seems rather well under- 
stood (see e.g. Pamyatnykh 1999). For low I modes theory 
predicts the occurrence of two instability domains within the 
main sequence band. In models of stars with earlier spectral 
types (Bl - B3) low order p- and g-modes are pulsation- 
ally unstable, as found in (5 Cep stars, and in stars with 
later spectral types (B4 - B8) high order g-modes are pre- 
dicted unstable, typically found in SPB stars, v Eri is per- 
haps the first j3 Cep star where both kinds of modes are 
present (Handler et al. 2004). In the following, we assume 
that the 0.432 cd _1 signal in the light curves of v Eri is an 
independent pulsation mode. There is a slight chance that 
it may be a combination signal, but its required order would 
be much higher than those of all the observed combination 
frequencies of similar amplitude. Thus we regard such an 
interpretation highly unlikely. 

In our fitted models we find unstable modes only in the 
4.0 — 6.5 cd _1 frequency range. The instability measure, rj, 
plotted in Fig. 7, is not a monotonic function of frequency. 
There is a noticeable maximum of the instability measure 
near v = 0.5 cd _1 for the I — 1 sequence and at somewhat 
higher frequency for 1 = 2. The existence of this maximum 
suggests that pulsational instability in this part of the eigen- 
mode spectrum may be found upon some modification of our 
models. 

Our proposal for extending the instability range is not 
original. We follow the idea of Charpinet et al. (1996, 1997) 
and Fontaine et al. (2003) that the combined effect of set- 
tling and radiative levitation may lead to significant over- 
abundance of iron in the driving zone. Since we do not have 
a stellar evolution code treating these effects, we simply in- 
troduced an ad hoc factor 4 enhancement of the iron-group 
elements in the log T range 5.1 — 5.5 and matched it smoothly 
to the standard solar value outside. The result is shown in 
Fig. 7. We may see that, while an extension of the instability 
to higher order p-modes is easy, obtaining instability in the 
high-order g-mode range is more difficult. We found unsta- 
ble i = 2 g-modes in the model calculated with a ov =0.1, 
which is cooler, but only at v « 0.8 cd _1 . Because of the 
ad hoc nature of the proposed iron enhancement we regard 
these results merely as an indication of a direction where the 
solution of the driving problem of modes detected in v Eri 
may be found. The proposal that the diffusion effects may 
play a role in this star seems not unreasonable in view of its 
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Figure 7. The instability measure, r/, for £ = 1 and 1 = 2 modes 
in a standard model calculated with a ov = and in the corre- 
sponding model with the Fc enhancement in the driving zone. 
rj > for unstable modes. 77 = — 1 if there is no active driving 
layer in the star for a particular mode and it is +1 if driving oc- 
curs in the whole interior. The position of the modes detected in 
v Eri are marked with bars on the rj = axis. 



very slow rotation which should not induce any appreciable 
mixing of elements and this is the first star where modes in 
such a broad frequency range were detected. However, slow 
rotation cannot be the sufficient condition for a broad-band 
excitation. In HD 129929, which is even a slower rotator, all 
modes detected so far are confined to the narrow frequency 
range of 6.45 - 7 cd _1 . 



5 POSSIBLE IDENTIFICATION OF THE 
REMAINING MODES 

Our manipulations on the Fe abundance in outer layers are 
not without consequences for mode frequencies and hence on 
the seismic model of v Eri. It is not difficult to find a model 
with modified Fe which fits the three frequencies used in 
section 2.1. The model cannot be very different; however, the 
frequencies of higher order p-modes are more significantly 
changed. 

Let us begin with the possible mode at v = 0.432 cd _1 . 
Though it is easier to find instability at higher degrees, the 
i = 1 identification seems more plausible because in this case 
the maximum of rj occurs closer to the observed frequency. 
Assuming i = 1 and the rotational splitting the same as in 
the (I = 1, g x ) mode, we obtained 0.429 and 0.435 cd _1 
for the frequency of the [I = l,m = — l,gie) mode in the 
models with a ov = and 0.1, respectively. Thus, we regard 
such an identification for the v = 0.432 cd" 1 peak as likely. 

For the observed variation at v = 7.20 cd~\ the 
only possible identification, regardless of the choice of a ov 
and modification in the Fe abundance, is a mode of the 
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(£ = 2,p ) quintuplet. In this case the effect of the Fe en- 
hancement is more significant. The enhancement described 
in the previous section leads to a frequency decrease by 0.02, 
which is comparable with the rotational splitting. With this 
model the best fit is for m = 1, while with the standard 
models it is for m = 2. 

The effect of the Fe enhancement increases with mode 
frequency. This is not surprising because the acoustic prop- 
agation zone expands to outer layers where the effect of 
the enhancement is largest. For the observed mode at v = 
7.90 cd~ we know the lvalue from photometry, and it is 1. 
In our standard seismic models the £ — l,p2 modes have fre- 
quencies 8.10 and 8.03 cd _1 in models with q ov = and 0.1, 
respectively. In this case, even with the freedom in choosing 
the m value we cannot reproduce the measured frequency. 
The frequency shift due to the Fe enhancement, which is 
about 0.1 cd _1 , brings the calculated frequencies close to 
the fit. We regard this fact as partial support for the pro- 
posed effect. 

We note that fitting the frequencies of all the modes 
detected in the 5.6 - 7.9 cd _1 frequency range is possible 
in standard evolutionary models if one allows lower effective 
temperature and overshooting in the range 0.2 < a ov < 0.3 
(Ausseloos et al. 2004). 



6 CONCLUSIONS AND DISCUSSION 

We believe that our seismic models of v Eridani yield a 
good approximation to its internal structure and rotation 
but there is room for improvement and a need for a full ex- 
planation of mode driving. We did not fully succeed in the 
interpretation of the oscillation spectrum of v Eridani with 
our standard evolutionary models and our linear nonadia- 
batic treatment of stellar oscillations. We fail to reproduce 
mode excitation in the broad frequency range, as observed, 
and to reproduce the frequency of the (I — 1,^2) mode, as- 
sociated with the highest frequency peak in the spectrum. 
We showed that both problems may be cured by allowing an 
enhancement of the iron abundance in the zone of the iron- 
opacity bump. The enhancement may result from effects of 
radiative levitation, as first proposed by Charpinet et al. 
(1996, 1997) to explain oscillations in sdB stars. Our pro- 
posal is based on models with an ad hoc iron enhancement 
in the bump zone and must be checked with use of stellar 
evolution codes that take into account effects of levitation 
and diffusion of chemical elements. We regard this as the 
most timely theoretical work aimed at the interpretation of 
v Eridani. 

The proposed modification concerns only the outer lay- 
ers and has a small effect on models of the stellar interior 
and the frequencies of g and pi modes. Therefore we believe 
that our seismic models describe the deep internal structure 
of the star. These models reproduce frequencies of the fun- 
damental radial modes and two 1=1 modes and have effec- 
tive temperatures and luminosities within the measurement 
error box. Satisfactory models exist in a range of the over- 
shooting parameter, q ov = (0 — 0.12), which corresponds to 
the mass range (9.9 — 9.0)Mq, age range (16 — 20) My, and a 
fractional hydrogen depletion in the core range (0.34 — 0.38). 
Here the important finding is the upper limit for the extent 
of overshooting. There is a prospect for getting a more strin- 



gent constraint on a ov with a more accurate determination 
of the star's effective temperature. 

We believe that our most important finding concerns 
internal rotation. We presented evidence for a sharp decline 
of the rotation rate, Q, through the /i-gradient zone, around 
the shrinking convective zone. With data on the splitting for 
only two I = 1 modes, our estimates had to rest on the as- 
sumed form of the fl(r) dependence. Since a significant gra- 
dient of £1 is possible only in the /i-gradient zone, we adopted 
a constant above it (rotation within the convective core 
has a negligible influence on the splitting) . With this simpli- 
fication we found that data require a decline by a factor ~ 5 
within the zone and an equatorial velocity of about 6 km/s. 
The second order effect from this slow rotation in the en- 
velope cannot account for the asymmetry of the [I = 1, g x ) 
triplet as determined from the pre-campaign data. Only with 
new observations we will be able to tell whether the prob- 
lem is real. Our star is not the first j3 Cephei star for which 
evidence for nonuniform rotation was put forward. However, 
we believe that in the v Eri case the evidence is significantly 
stronger than in the case of HD 129929 (Aerts et al., 2003), 
because we relied on splitting data for modes having very 
different probing kernels. 

v Eridani proved to be a very important pulsating star. 
After the Sun, it so far is perhaps the most rewarding main 
sequence object for asteroseismology. Many more modes 
were detected in several S Scuti stars but in none of them so 
many modes have unambiguously been identified. This star 
deserves continued observational efforts. An extended data 
base is needed for precise determination of the frequencies 
of the modes in the 1 = 1 triplets. This is important for a 
more credible assessment of differential rotation. In this con- 
text it would be important to obtain a precise value of the 
projected equatorial velocity of rotation from spectroscopy. 
The best v sin i value currently available (Aerts et al. 2004) is 
higher than that of equatorial velocity inferred in this work 
which may however be in part due to the combined effects 
of pulsational and thermal line broadening. 
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